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We discuss the applicability of the Boltzmann equation to the classical two-dimensional electron 
gas. We show that in the presence of both the electron-impurity and electron-electron scattering 
the Boltzmann equation can be inapplicable and the correct result for conductivity can be different 
from the one obtained from the kinetic equation by a logarithmically large factor. 



The description of kinetics of the classical electron gas 
in semiconductors is usually done with the help of the 
Boltzmann kinetic equation Jy 



dn p (r,t) dn p (r,t) 



eE 
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dt dr dp 

Iimp(n> P (r,t)) + I ee (n p (r,t)) + I p h(n p (r,t)) 
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where n p (r,t) is the electron distribution function, E 
is the external electric field, v and p are the electron 
velocity and momentum, and limp, lee and I p h are the 
electron-impurity, electron-electron and electron-phonon 
scattering integrals which are characterized by the relax- 
ation times Ti, r e and r p h respectively. 

At low enough temperatures T when Tj <C T p h, T e , the 
conductivity of the sample is determined by the electron- 
impurity scatterings and is given by the Drude formula 



(Jq = e Dv 



(2) 



Here D = vpTi/2 (where vp is the Fermi velocity, 
Ti = l/upA^a, Ni is the impurity concentration and a is 
the amplitude of the electron-impurity scattering), and 
v is the electron density of states at the Fermi level. To 
be concrete, we consider the case of a degenerate electron 
gas, though our results are of a general character. Eq. (|^) 
also holds in the case of intermediate temperatures when 
Tph 3> Ti 3> T e . Indeed, since the electron-electron scat- 
terings conserve the total momentum in the electron sys- 
tem, they do not contribute directly to the resistance of 
the system. 

The form of limp in Eq. (|l]) is a result of a procedure in 
which averaging over realizations of the random impurity 
potential is done before solving the equation of motion. 
The criteria of applicability of the conventional Boltz- 
mann kinetic equation (|l|) correspond to the absence of 
correlations between subsequent scattering events. This 
leads us to a requirement that the electron gas is weakly 
interacting 



Ti,T e > t c 



(3) 



where t co u is the the collision time, which is estimated 
as a/v in classical mechanics, and as | in quantum me- 
chanics (e is the characteristic electron energy). 



It has been shown g that in the two-dimensional case 
even under conditions the quantum interference cor- 
rections to the Drude formula logarithmically diverge as 
temperature T goes to zero. This means that in this case 
the kinetic equation (|l|) is not applicable at small T. 

In this article we study the case of high enough tem- 
peratures when the quantum interference corrections can 
be neglected. The classical corrections to Eq. (||) asso- 
ciated with correlations between subsequent scatterings 
were considered in in the framework of the Lorentz 
gas model. It was shown that in the two-dimensional case 
correlational corrections to the diffusion coefficient and 
to the conductivity are small, but nonanalytical functions 
of the concentration of scatterers iVj. 

We consider in this article the case of a two- 
dimensional electron gas with the mean free path l e — 
VFT e much smaller than the average distance between 

— 1/2 

the impurities: l e <C N t . We show that under condi- 
tions (|3|) a combination of electron-electron and electron- 
impurity scatterings can lead to big correlational correc- 
tions to the Drude formula and that the correct expres- 
sion for the conductivity of the system is 



a = (T (l 
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where r\ is the kinematic viscosity of the electron liquid, 
which is estimated for a Fermi liquid as r\ ~ vpl e . Eq. (^) 
differs from Eq. (^) by a logarithmically large factor, pro- 
vided 
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> 1. 



On distances from impurities larger than l e , the motion 
of the electron system can be described by the Navier- 
Stokes hydrodynamic equation which can be derived from 
Eq. ([[]) with limp = H . In this case the elastic scatter- 
ers should be described by a random potential and the 
averaging over the random realizations of the potential 
should be done after solving the Navier-Stokes equation. 
This should in principle give a result for conductivity that 
is different from the one implied by Eq. ([!]). The two 
approaches are equivalent in the approximation which 
assumes the variation of the hydrodynamic velocity be- 
tween impurities to be small. 



f 



The idea of this paper is that, in fact, in the two- 
dimensional case the hydrodynamic velocity of electrons 
is a significantly nonuniform function of coordinates: the 
velocity near impurities is much smaller than inbetween 
them. This is in contrast to the three-dimensional case 
where the spatial variation of the velocity is small and 
Eqs. (0), (||) are applicable. The reason for the signifi- 
cant nonuniformity of the velocity field is connected to 
the Stokes paradox [|| : in the presence of an obstacle by- 
passed by the fluid, the solution of the linearized Navier- 
Stokes equation in the two-dimensional case grows log- 
arithmically with the distance from the obstacle. This 
means the absence of a bounded solution. In our case 
the finite concentration of scatterers provides a bound on 
the magnitude of the velocity field. However, the magni- 
tude of the average velocity of electrons induced by the 
electric field turns out to be logarithmically bigger than 
that which follows from Eq. ([!]). 

In the approximation linear in E, the Navier-Stokes 
equation has the form 



Vp + uAu + neE - ^ fS(r - n 
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Here n is the electron concentration, \x = mnr\ is the dy- 
namic viscosity of the electron fluid (where m is the elec- 
tron mass), and r$ are coordinates of impurities labeled 
by the index "i". The last term in Eq. (^) is associated 
with the force on the fluid by the i-th impurity, whose ra- 
dius is much smaller than l e . Therefore the value of the 
force / can be estimated assuming that to zeroth order 
in a/l e the distribution function of the electrons near the 
i-th impurity (|r— Ti\ < l e ) is an equilibrium one with the 
average velocity m(t*j). (In the following we will assume 
that all u(ri) ~ Uo are of the same order.) Consequently 
we have l| 



/ = np F au . 



(6) 



In the stationary case the conservation of the total mo- 
mentum of the system gives 



eE = fNi/n. 



(7) 



Far from the position of impurities (\r — S> l e ) the 
electron liquid can be considered as an incompressible 
one. The continuity equation for an incompressible liq- 
uid 

V u = 

and the Navier-Stokes equation (||) have a solution for 
the pressure and the velocity 



u(r) = Cx + x— lnr — — (f cos( 



P(r) =Po- MCo 



cos 9 



(8) 



where Cn = tJ—. f and 9 are unit basis vectors of the 

u Z7TfL 1 

cylindrical coordinate system and the x-axis (with the 
unit vector x) is oriented parallel to the electric field. 

The value of the coefficient C is determined using 
u(r = l e , 9 = tt) = M and Eqs. (|), (§) and (§): 

<J E av F 

C = (1 - - — lnZ e ). 
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As a result, the spatially averaged velocity u of the 
electronic flow is 
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n e 
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Expressing the current density as j = neu, the electric 
conductivity Eq. (^) is obtained. 

The calculations presented above assumed that the vis- 
cosity of the electron gas is a well-defined quantity. On 
the other hand it is well known that the fluctuational 
corrections to the viscosity of a 2D neutral liquid diverge 
0, which means that hydrodynamic equations are non- 
local in two dimensions. This fact was pointed out in j^] 
and discussed in the most complete form in Below 
we present a trivial generalization of the procedure [^) to 
estimate fluctuational corrections to the viscosity of the 
two-dimensional charged electron liquid. In this case the 
low-lying collective excitations are plasmons and the fluc- 
tuations of the charge and of the hydrodynamic velocity 
can be described as a combination of plasmon excitations 
with different wave vectors q. The real par t of the plas- 
mon frequency is |Io[ Lj q = ^J c lnne 1 \q\jm , where q is 
the plasmon momentum. Since the imaginary part of the 
frequency uo' q = i]\q\ 2 /2 is proportional to \q\ 2 (77 = [ijmn 
is the kinematic viscosity of the liquid) , fluctuations with 
sufficiently low \q\ have arbitrarily large mean free path 
and therefore plasmons are well-defined excitations. This 
is the physical reason why the fluctuational mechanism 
of viscosity is always the basic one at sufficiently small 
gradients and in samples with no elastically scattering 
potential. Since the phases of the thermal fluctuations 
are random it is convenient to introduce a distribution 
function of plasmon fluctuations n(q), normalized in such 
a way that the quantity J j0y I u(q)n(q) is the average 
energy density of fluctuations. The flux of momentum 
associated with these excitations is [[fl| 



Uij = - 



dq duj„ 



(9) 



The part of the distribution function of plasmons which 
is proportional to the gradient of the hydrodynamic ve- 
locity can be retrieved from the kinetic equation 



dm q dn u (q) 
dq 



Or 



-(n(r,q) - n u {q))/r 



(10) 



w * is the plasmon mean free time and n u (q) 
is the equilibrium Bose distribution function correspond- 
ing to a local hydrodynamic velocity u(r). At small ui 



2 



and \u\ one can approximate n u (q) ~ - _g. u ^ r y Using 
Eq. ( |Tc| ) and expanding n u (q) in it, one can estimate 
from Eq. (|) the fluctuational contribution to the viscos- 
ity 



1 f dui q 2 <9n (<?) 
ran / Off ouj a 



(11) 



(where rto(<z) == n u= o(q)). Thus the fluctuational cor- 
rection to the viscosity can be described only in terms 
of hydrodynamic quantities: viscosity and plasmon dis- 
persion relation. In two-dimensional pure systems the 
expression Eq. ( pi] ) logarithmically diverges. However, 
in the presence of the elastically scattering potential the 
dispersion relation of plasmons which we used holds only 
at those frequencies for which the plasmon mean free 
time is shorter than the electron-impurity scattering time 
t < Ti. This determines a cut-off in Eq. (p"l|). As a re- 
sult, the fluctuational correction to the viscosity can be 
neglected provided 

T T (To . 
6r) m(- jtt) < V- 

In conclusion, we would like to stress that the exis- 
tence of the large logarithmic factor in the Eq. (^) and 
the consequent inapplicability of the Boltzmann kinetic 
equation is a property of the two-dimensional problem. 
In the three-dimensional case the Stokes paradox does 
not exist, the difference between Uo and u is small and 
the Boltzmann kinetic equation can be applied. 
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